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Abstract. Periodic light curve behavior predicted for some binary black hole sys- 
tems might be detected in large samples, such as the multi-million quasar sample ex- 
pected from the Large Synoptic Survey Telescope (LSST). We investigate the false- 
alarm probability for the discovery of a periodic signal in light curves simulated using 
damped random walk (DRW) model. This model provides a good description of ob- 
served light curves, and does not include periodic behavior. We used the Lomb-Scargle 
periodogram to search for a periodic signal in a million simulated light curves that 
Q . properly sample the DRW parameter space, and the LSST cadence space. We find that 

?— i ■ even a very conservative threshold for the false-alarm probability still yields thousands 

c/3 | of "good" binary black hole candidates. We conclude that the future claims for bi- 

nary black holes based on Lomb-Scargle analysis of LSST light curves will have to be 
interpreted with caution. 
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1. Introduction 

Modern surveys of the sky, such as the Sloan Digital Sky Survey (SDSS, lYork et all 
2000h . have collected huge amounts of data (20 TB for SDSS), requiring the develop- 
ment of automate d analysis methods. The Large Synoptic Survey Telescope (LSST, 



Ivezic et al.ll2008l) will gather even more data than the SDSS did (one SDSS equivalent 
per night over ten years of operations). Among other populations, LSST will identify 
several million quasars (QSO) and obtain their light curves. In this contribution, we 
discuss an automated analysis of a million simulated light curves to search for periodic 
variability. 

The optical variability o f QSOs has been recognized since they were first identified 
(Matt hews & Sandagd fl963) and is aperiodi c and on the order of 2 0% on timescales of 
months to years (for recent results see, e.g., iMacLeod et aLll2010l hereafter M10). Pe- 
riodic variability has bee n suggested as one of the observational char acteristics of a 
binary black hole system (lKomossall2003l) . (also see lShen & Loebll2010l. for an investi- 
gation of broad line emission spectra in binary black holes), but there is no convincing 
observational evidence for such systems yet. 
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2. Goal 



The large QSO sample expected from LSST might enable a detection of periodic sig- 
nal in observ ed light curves . Recently we employed the Lomb-Scargle periodogram 
(lLomblll976l : IScarglel[l982l : iHorne & BaliunaslTT986h to test ~9000 spectroscopically 
confirmed QSOs from SDSS Stripe 82 (S82) for periodic variability (see Appendix of 
M10). We reuse the tools developed for the analysis of SDSS data and study light 
curves simulated using a mathematical variability model trained on SDSS data, and 
sampled using simulated LSST cadences. 



3. Creation of Light Curves 

Approximately lxlO 6 QSO light curves were gener ated using the damped random 
walk model (Kell y et al.ll2009l : iKozlowski et al.1 20101. M10). The difference between 



this model and the well-known random walk is that an additional self-correcting term 
pushes any deviations back towards the mean flux on a time scale t. The above studies 
have established that DRW can statistically explain the observed light curves of quasars 
at an impressive fidelity level (0.01-0.02 mag). 

The input parameters to the model are the characteristic time scale, r, and the root- 
mean-square (rms) variability on long time scales, or structure function, sf. The input 
parameters were determined using the scalings with black hole mass (Mbh), absolute 
magnitude (Mi), and redshift found by M10. These physical parameters were drawn 
from the distribution shown in Figure 1. After generating ~83,500 well-sampled light 
curves, each light curve wa s res ampled to the 12 different simulated r-band LSST ca- 
dences from Delg ado et all (120061 ~200 observations spread over 10 years) to obtain 
~ lxlO 6 total light curves. 

All light curves were analyzed with Lomb-Scargle periodogram software. If the 
maximum power spectral density (PSD) value was above the level set according to 
the false alarm probability, fap, then the curve became part of a set used for further 
examination. 



4. Results 

We found that for high fap values, the actual number of light curves which exceeded 
the fap level was less than the theoretical value by an appreciable amount. For a fap 
value of 5%, there were 13,294 (1.3%) light curves that exceeded that level, for afap 
value of 1%, there were 4696 (0.47%), and for a fap value of 0.1%, there were 1035 
(0.1%). 

Figure 2 shows that values of r for QSO models with PSD exceeding the 1% fap 
level (red histogram) are distributed differently than those for all the QSO models (black 
histogram). This bias is due to the fact that when t is long, only a few "oscillations" 
are observed over the duration of the light curve, causing the periodogram to mistake 
the damped random walk behavior for a periodic behavior. 

Figure 3 compares the proportions of QSO models found to exceed the 1% fap 
level for each of the cadences. There is a marked difference between the various ca- 
dences here, showing that some of the cadences allow many more periodogram results 
for which the \%fap level is exceeded. This may be because the variations in the sam- 
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Figure 1 . QSO model parameter ranges used to derive inputs to light-curve generation 
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Figure 2. Distribution of characteristic time, r, for all QSO light curves compared 
with t distribution for light curves with a peak SED value exceeding the \ %fap level 



pling pattern and the "windowing" effect of the Lomb-Scargle periodogram can cause 
a resonance effect at some test frequencies. Therefore, the results obtained by using 
the Lomb-Scargle periodogram method can be greatly influenced by the pattern of the 
observations used. 



4 



D. Westman, C. MacLeod, Z. Ivezic 




% of total (4696) OSOs for each cadence □ 1 2 



Figure 3. Proportion of QSOs with peak PSD values over \ %fap level sorted by 
LSST cadence 



5. Conclusion 

This work shows that the Lomb-Scargle periodogram method may be useful for detect- 
ing potentially periodic behavior in QSO light curves in a large-scale surveys, such as 
the one to be carried out by the LSST. However, even with fap as small as 0.1%, the 
large LSST sample would yield ~ 1,000 false candidates. Therefore, future claims of 
periodic behavior based on Lomb-Scargle analysis of LSST light curves will have to be 
interpreted with caution. In particular, black hole binary candidates identified by this 
method would have to be examined individually and with supplemental observations. 

Acknowledgments. We acknowledge support by NSF grant AST-0807500 to the 
University of Washington, and NSF grant AST-0551161 to the LSST for design and 
development activity. Z. Ivezic thanks the University of Zagreb, where a portion of 
this work was completed, for its hospitality, and acknowledges support by the Croatian 
National Science Foundation grant 0-1 548-2009. 

References 

Delgado, E, et al. 2006, in Observatory Operations: Strategies, Processes, and Systems, edited 
by D. R. Silva, & R. E. Doxsey (Bellingham, WA, USA: SPIE), vol. 6270 of Proc. SPIE, 
6270 ID 

Home, J. H, & Baliunas, S. L. 1986, ApJ, 302, 757 

Ivezic, Z., Tyson, J. A., Allsman, R., Andrew, J., Angel, R., & for the LSST Collaboration 

2008, ArXiv e-prints. larXiv: 0805 ■ 23661 
Kelly, B. C, Bechtold, J., & Siemiginowska, A. 2009, ApJ, 698, 895 

Komossa, S. 2003, in The Astrophysics of Gravitational Wave Sources, edited by J. Centrella 

(Mellville, NY, USA: AIP), vol. 686 of AIP Conf. Proc, 161 
Kozlowski, S., et al. 2010, ApJ, 708, 927 
Lomb, N. R. 1976, Ap&SS, 39, 447 
MacLeod, C. L., et al. 2010, ApJ, 721, 1014 
Matthews, T. A., & Sandage, A. R. 1963, ApJ, 138, 30 
Scargle, J. D. 1982, ApJ, 263, 835 

Shen, Y, & Loeb, A. 2010, ApJ, 725, 249 

York, D. G., et al. 2000, AJ, 120, 1579. |arXiv:astro-ph /iS)<S) 6)6396 



